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Abstract 

Deviations from randomicity in In x Gai_ x As semiconductor alloys induced 
by elastic effects are investigated within the Keating potential. Our model is 
based on Monte Carlo simulations on large (4096 atoms) supercells, performed 
with two types of boundary conditions: Fully periodic boundary conditions 
represent the bulk, while periodic boundary conditions along the x and y 
directions and a free surface in the z direction simulate the epitaxial growth 
environment. We show that In-In correlations identified in the bulk tend to 
be enhanced in the epitaxially grown samples. 
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I. INTRODUCTION 



In recent years, order versus disorder in substitutional semiconductor alloys has moti- 
vated several theoretically and experiment alEHLll studies. Special emphasis is given to the 
understanding of the physical mechanisms behind observed deviations from random distri- 
butions in the atomic positions. The distribution of different species over the atomic sites 
in such compounds is largely responsible for variations in electronic properties such as the 
band gap, the density of states at the Fermi energy and electron confinement levels. More- 
over, the atomic correlations in these systems are responsible for processes such as interface 
segregation and clustering, which are fundamental in defining the atomic scale structure and 
roughness of heterostructure interfaces. 

In particular, In^Gai-aAs semiconductor alloys with x < 20% have been the sub- 
ject of experimental studies indicating deviations from randomicity in their atomic 
configurations!^!. Zheng et a/.i identified a strong correlation between In atoms in x = 20% 
samples. Clusters of 2 - 3 In atoms were reported to form preferentially along the [001] 
growth direction, in second-nearest-neighbor (2nn) positions. Chao et a/.i investigated 
x = 5% samples and, in apparent contradiction with Zheng's results, no strong correlations 
were found in for 2nn pairs along the growth direction, whereas important anti-correlation 
for first-nearest-neighbor (Inn) pairs along [110] was reportedEl. The absence of clustering 
of In atoms in x = 12% quantum wires was also noted by Pfister et a/.i . Chao and Zheng 
emphasize that the mechanism behind the clustering/anticlustering effects in these alloys 
is the strain induced by In atoms, due to the 7% mismatch between the GaAs and InAs 
lattice parameters. Both studiesii were based on In^Gai-^As alloys grown using molecular 
beam epitaxy (MBE) on GaAs (001) substrates at temperatures ~ 800f^. The alloy samples 
were kept below the critical thickness on GaAs substrate. Cross-sectional images for cleaved 
(110)0 and (lio! surfaces, obtained from scanning tunneling microscopy (STM) were used 
to determine the In distribution in the alloy. We follow here the orientation convention of 
Pashley et a/H. 
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The only apparent reason for the discrepancy concerning the second-neighbor clustering 
tendency in these experiments is the In concentration range: Cluster formation was reported 
m x = 20% sampleJ, while no second-neighbor correlation was identified in samples with 
smaller values of x&i. We investigate this possibility presenting a theoretical study for In- In 
pair correlations within a Keating model potential^ to describe the elastic interactions in 
the alloy. The Monte Carlo Metropolis algorithm is used to determine the thermodynamic 
equilibrium configuration at finite temperatures. Phase-diagram calculations predict that 
InGaAs alloys are completely miscible at the growth temperature!!, but no detailed theo- 
retical analysis of pair correlation functions has been found in the literature. Therefore, 
one of the goals of the present work is to determine if the resulting In-In correlations are 
a consequence of bulk thermodynamics or if growth kinetics is important. Our simulations 
are then performed within two models, according to the boundary conditions: A bulk model 
(fully periodic) and an epitaxial growth model (periodic only in the x and y directions). 



II. FORMALISM 



Following Ref.i, we define the pair correlation function, C (r 12 ) as 
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where -P(r 12 ) = P(r 1; r 2 ) is the probability of finding two In atoms at and r 2 , which, 
for homogeneous alloys, is only a function of the relative position r 12 = r x — r 2 , and i?(r 12 ) 
is the equivalent quantity calculated for a random alloy with the same composition. If 
C (i"i 2 ) > 1, we have correlation between the atoms that form the pair, and if C (i"i 2 ) < 1, 
we have anticorrelation. Correlation indicates an effective attraction between the atoms, 
and anticorrelation implies a net repulsion. 

The Keating potential provides a good description for elastic energies in III-V semicon- 
ductor alloysEl'0. It is given by 



r 2 - - d-- 

1 13 U V 



E— — 

lijk 8di i d i k 



^3 



r jk + 



dijdjh 



(2) 



where dij is the equilibrium bond length between atoms i and j in the corresponding 
zincblend binary compound, r^- is the relative position vector between nearest neighbors in 
different sublattices and o;^ and (3^ are the bond-stretching and bond-bending constants, 
respectively. The elastic energy thus results from two contributions: The first summation 
(performed over nearest-neighbor pairs) refers to the excess energy due to the bond-length 
variations and the second term (three-center term) gives the contribution due to the devia- 
tions of bond angles from their ideal tetrahedral values. Values for d^, otij and flijk were taken 
from Ref.0. In particular, the difference between <icaAs = 2.448 A and di n A s = 2.622 A leads 
to the structural mismatch in the system. For heterogeneous bond angles, flijk was taken as 
the geometric mean of the bond-bending constants for the binary compounds: 



Calculations are carried out in supercells with N x N y N z conventional cubic unit cells 
of the fee lattice along the x, y and z directions, respectively, resulting in a system with 
iV = 8N x N y N z atoms. In^Ga^^As semiconductor alloys are modeled by supercells in which 
x ■ N/2 sites of the group-Ill sublattice are occupied by In atoms and the remaining sites 
by Ga, while the N/2 group-V sublattice sites are occupied by As atoms. For a given 
configuration, the energy is minimized by relaxing all the atomic positions in the supercell. 
Initially, we set the atoms to be randomly distributed, and evolution to thermodynamic 
equilibrium at the growth temperature (T = 800 K) is carried out through the Monte Carlo 
Metropolis algorithm, via first-neighbor In <-> Ga positions exchanges. The calculations are 
performed at constant volume. In the bulk model (Section 3) we use the Vegard's law! to 
give the approximated value for the lattice parameter and for the epitaxial growth model 
(Section 4) we keep the lattice parameter to the GaAs value, aiming to describe a strained 
growth process. 




(3) 
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III. BULK MODEL 



For the bulk model, after ~ 10 3 MC steps the system attains thermodynamic equilibrium 
and statistically independent (~ 500 MC steps apart) alloy configurations are collected into 
an ensemble to calculate the pair correlation function. This procedure is repeated for 6 initial 
random alloy configurations, leading to analyzed ensembles of at least 100 configurations for 
each temperature and composition. 

The third column of Table | shows our results for x = 0.05 and x = 0.20 in In x Gai_ x As. 
Notice an anticorrelation (repulsion) between Inn pairs and a correlation (attraction) be- 
tween 2nn pairs. These results are in qualitative agreement with experiments. However, 
quantitative agreement is poor: Chao'si experiments indicate a much stronger anticorrela- 
tion for Inn pairs than what is predicted by bulk thermodynamics. In addition, the observed 
experimental discrepancy regarding the 2nn correlation cannot be explained as a composi- 
tion effect within the bulk model: There is no increase in the calculated correlation upon 
increasing In concentration from 5% to 20%. 

Although not significantly larger than one, the 2nn correlation is clearly present. It is 
somewhat surprising that there can be such an effective attraction (for 2nn pairs) between 
In atoms based on elastic effects alone. Since In atoms are larger than Ga ones, one should 
expect that In impurities should compress the lattice locally and therefore repel each other. 
To see this in more detail, let's consider the interaction of two isolated In impurities in GaAs. 
We define a pair interaction energy as 

AE = U- 2E , (4) 

where U is the Keating energy for the supercell with two impurities and Eq the elastic energy 
for a single impurity, calculated to be Eq = 138 meV. The interaction energy is displayed 
in Table ||. Notice a net attraction (negative energy) along the [001] direction and a net 
repulsion (positive energy) along [110], consistent with the MC results. The interaction 
energy is therefore highly anisotropic, as can also be visualized in Figure |1| These results 
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can in fact be explained by the intricate geometry of the zincblende structure, as we can see 
by considering the mean Ga-As bond length deviations induced at nearby site by a single 
In impurity 

1 4 

ALj = - \ r ij \ ~ rf GaAs , (5) 

4 i=i 

where % is a Ga site and the sum in j is performed over the As sites nearest neighbors 
to i. Negative (positive) AL at a given site means that the lattice is locally compressed 
(expanded) at that site by the nearby In impurity and therefore suggests that putting a 
second In impurity at that site will be energetically unfavorable (favorable). One can see 
from Table |I| that even a larger impurity such as In in GaAs can produce a local expansion of 
the zincblende structure along certain directions (for instance, the [001] direction). An exact 
opposition between the signs of AL and AE can also be readily seen ion the Table, thus 
supporting the argument given above. The highly anisotropic nature of the pair interactions 
is in agreement with first-principle calculations for these energies^ Interaction energies are 
determined there from a series expansion in terms of multiatom contributions, fitted to a set 
of ordered structures. Their results are in qualitative agreement with the present Keating 
model calculations, giving the most repulsive interactions between first and fourth neighbors 
[(110) and (220)], while the less repulsive ones correspond to second and third neighbors 
[(002) and (211) respectively]. Quantitatively, our results in Table [II] are lower than those 
in Ref.i by several meV and, in particular, no attractive pair interactions result from the 
first-principle fits. 

IV. THE EPITAXIAL GROWTH MODEL 

Aiming at a more realistic model describing the experiments, we propose a simple de- 
scription of epitaxial growth: (i) No periodic boundary conditions are imposed in the z 
(growth) direction; (ii) The first two (bottom) monolayers, with N x = N y = 8, represent 
the GaAs substrate, and the atomic species are fixed as in bulk GaAs. Atomic positions are 
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fixed according to GaAs only in the z = plane; (in) Growth in the z direction proceeds 
by adding a single alloy monolayer at a time, until N z = 8: The In <-> Ga exchanges in 
the MC algorithm are restricted to the most recently grown monolayer (free surface) and 
the next monolayer is introduced only after thermal equilibrium has been attained. Atomic 
correlations between successive monolayers are thus induced by the In distribution in the 
inner layers. This is a reasonable assumption since, at typical MBE growth rates, the time 
scales for atomic processes at the free surface are much shorter than those in the previously 



grown planes. In addition, as can be seen in Table |HI|, there is an elastic-induced attraction 
between In impurities and the surface. Allowing bulk-surface In exchanges would produce 
In segregation at the growing surface, an effect that is known to occur but which we do not 
intend to describe here. In this present model, we do not treat dangling bonds, i.e., the 
surface is described simply as a truncated semi-infinite solid. 

In the fourth column of Table p], we present the pair correlation function for 5% and 
20% In concentrations in In^Gai-^As. There are now two inequivalent in-plane Inn pairs 
directions, [110] and [110]. For the epitaxial growth model, the configurational averages 
were performed over 20 different configurations for each concentration and temperature. 
The statistical ensemble in the epitaxial growth model is considerably smaller than for 
the bulk model, since each growth process here contributes with a single configuration to 
the ensemble. Moreover, the bottom monolayers are fixed in composition to be of GaAs, 
which decreases the total number of In pairs in each generated configuration. This leads 
correspondingly to larger statistical error bars. 

We note that for x = 5%, C (r 12 ) along [110] is significantly reduced with respect to the 
bulk model behavior (see Table |), while pairs along [110] become essentially uncorrelated. 
This is an interesting feature of the epitaxial growth model: The free surface induces impor- 
tant anisotropy between the [110] and [110] directions, and the correlation function in the 
bulk is approximately the average between those calculated for the two inequivalent direc- 
tions in the epitaxial growth model. The experimental results of Ref.i give C(t 12) = 0.25 
along [110] for x = 5% and T = 800 K. Our corresponding values for the bulk and epitaxial 



models, C(r 12 ) = 0.83 and 0.56 respectively, indicate that epitaxial growth is a decisive 
ingredient in the observed experimental behavior. Although the agreement is still not quan- 
titative, our results confirm that the observed Inn anti-correlation has an elastic origin. For 
2nn pairs with x = 20%, we do not observe significant increase of the correlation function 
above the random value. However, the correlation enhancement in going from 5% to 20% 
seems to be present in the epitaxial growth model, although error bars are too large to allow 
a more precise statement. Therefore, only elastic effects and a simple epitaxial growth do 
not provide a clear and quantitative explanation for the clusters formation reported in Ref.@, 
even though we detect a tendency for In-In attraction in the growth direction in many of the 
cases studied. Perhaps a more sophisticated growth model, including aspects such as surface 
reconstruction and step-mediated growth, is necessary to better evaluate such results. 
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TABLES 



TABLE I. 



Calculated (bulk and epitaxial growth models) and experimental correlation functions cal- 
culated for different In concentrations (x) in the alloy and In-In relative positions (r) in an 
undistorted GaAs supercell in units of half of the conventional GaAs lattice parameter a;. 



X 


r (x|) 


Bulk 


Epitaxial growth 


Experiment 


0.05 


(110) 


0.83 ± 0.03 


0.56 ±0.11 


0.251 




(110) 


0.83 ± 0.03 


0.93 ±0.09 






(002) 


1.10 ±0.05 


0.89 ±0.12 




0.20 


(110) 


0.84 ±0.01 


0.86 ±0.03 






(110) 


0.84 ±0.01 


0.96 ±0.02 






(002) 


1.07 ±0.01 


1.03 ±0.03 


>10 
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TABLE II. 



Elastic properties of the In-In impurity pair in GaAs. The columns give, respectively: (i) 
The relative position In-In positions (see Table I); (ii) The interaction energy; (iii) Bond 
length deviations induced by one of the impurities. 



r(x§) 


AE (meV) 


AL (10~ 3 A) 


(110) 


6.194 


-2.776 


(220) 


6.665 


-1.407 


(330) 


1.952 


-0.430 


(002) 


-4.077 


0.713 


(004) 


-0.637 


0.072 


(211) 


-2.869 


0.468 


(222) 


-0.939 


0.128 


(321) 


0.447 


-0.140 







TABLE III. 

Single In impurity energy as a function of depth for a free-surface supercell. Here, n indicates 
the monolayer where the impurity is: n = is the free (top) monolayer, n — 1 is the first 
monolayer below it and so on. The impurity energy increases with the distance from the 
free surface, meaning an elastic-induced effective attraction between impurity and surface. 
The energy converges to the bulk value at the third monolayer below the surface. 



n 





1 


2 


3 


Energy (meV) 


22 


105 


135 


138 
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FIG. 1. In-In elastic energy distribution in the (110) plane. Dark (light) regions correspond to 
lower (higher) energies. Open circles indicate the atomic sites in the group-Ill sublattice. In this 
plane, the [001] direction favors In pair formation and the [110] direction does not. 
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